Synergy in spreading processes: from exploitative to explorative foraging strategies 
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An epidemiological model which incorporates synergistic effects that allow the infectivity and/or 
susceptibility of hosts to be dependent on the number of infected neighbours is proposed. Construc- 
tive synergy induces an exploitative behaviour which results in a rapid invasion that infects a large 
number of hosts. Interfering synergy leads to a slower and sparser explorative foraging strategy 
that traverses larger distances by infecting fewer hosts. The model can be mapped to a dynamical 
bond-percolation with spatial correlations that affect the mechanism of spread but do not influence 
the critical behaviour of epidemics. 
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The identification of criteria to predict whether or 
not a spreading agent such as an infectious pathogen, 
rumour or opinion will invade a population is of great 
relevance both in biology and social sciences^"— li^. 
Numerous models have been proposed to gauge the 
invasiveness of spreading agents and assess the ef- 
fectiveness of control in preventing invasion of epi- 
dcmicsii^i2ii£iHiii"i2i. Within the context of infectious 
diseases, much theoretical work has been done for well- 
mixed populations but invasion in stochastic, spatially- 
structured, individual-based models^"— closer to real 
epidemics^^"— has also been considered. These mod- 
els do not deal with synergistic effects in transmission 
of infection thus assuming independent and identical 
action between hosts. This would suggest that mul- 
tiple challenges to a susceptible host from one, two 
or more neighbouring infected hosts are independent 
and not influenced by the local environment. How- 
ever, there is evidence for the existence of such effects 
in systems subject to colonisation by fungal and bac- 
terial pathogensiSiii, as well as in tumour growth2,iH. 
Synergistic effects have also been experimentally re- 
ported in studies of opinion dynamics^^, spread of be- 
haviour—, and animal invasion2i21. The model pre- 
sented in Ref. [l^ incorporates some temporal synergis- 
tic effects but it deals with well-mixed populations and 
spatial synergistic effects are not considered. Models 
for opinion dynamics and animal invasion have consid- 
ered some constructive synergistic spatial effects (e.g. 
population pressure^ and social impac t^^'^^ ). However, 
these effects are too simple to capture, for instance, 
possible changes in the foraging strategies of spreading 
agents that can significantly affect important features 
of invasions such as their size and time scales. Here, 
we present a model for spread of infection in spatially- 
structured populations and show that synergistic ef- 
fects in transmission of infection have non-trivial and 
significant consequences on epidemics. 

We consider an epidemic spreading through a pop- 
ulation of susceptible hosts placed on the sites/nodes 
of a square lattice of size L x L. The infection trans- 
mission rate between any donor-recipient (d-r) pair of 
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FIG. 1. Mechanism for (a) r-synergy and (b) d-synergy. 
Dark (red online) and lighter (green online) circles corre- 
spond to infected and susceptible hosts, respectively. Ar- 
rows indicate possible events for transmission of infection. 
For r-synergy, the susceptibility of a recipient susceptible 
host @ is enhanced (constructive synergy) or diminished 
(interfering synergy) in response to multiple simultaneous 
challenges from two or more neighbouring infected hosts 
(donors). The susceptibility of (?) depends on the number 
rir of neighbours simultaneously challenging it (wr = 3 in 
this example) . For d-synergy, the infectivity of a donor host 
(d) depends upon the number of connections of this host 
to other infectious hosts that can share resources with the 
donor host. In this example, rid = 1 (connection indicated 
by an edge). 



hosts depends on the number of infected hosts in the 
neighbourhood of the d-r pair. We focus on two par- 
ticular cases of the model denoted as r-synergy and 
d-synergy (Fig. [T]). 

The model is an extension of a basic spatial model 
for the SIR epidemic process in which hosts (sites) can 
be in one of the three states^: susceptible (S), infected 
(I) or removed and fully immune to further infection 
(R). Once a host is infected, it stays in such a state for 
a fixed unit of time, r = 1 and can pass infection dur- 
ing this infectious period to other S-neighbours, and 
then it is removed/recovered (I — >■ R transition). The 
infection process (S — !• I transition) occurs through the 
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pathogen being transmitted randomly with rate Ad-r(i) 
from one of the I-neighbours (donor) of the S-host (re- 
cipient). Synergistic effects make the rate Xd-r{t) de- 
pendent on the state of hosts in the neighbourhood of 
the particular d-r pair. By definition, Ad-r(i), is zero 
before the time of infection of the donor. Once the 
donor is infected, Ad-r(i) becomes positive but, in con- 
trast to the basic SIR modcl^ii, it can vary during the 
whole infectious period r due to possible changes in the 
neighbourhood of the d-r pair. The transmission rate 
can be conveniently split into two contributions, 

Ad-r(0 = max{0,a /3d-r(0} , (1) 

where a > is the elementary rate of infection for 
an isolated d-r pair exhibiting no synergistic effects. It 
does not vary over the infectious period and is assumed 
to be the same for all d-r pairs. The rate /3d-r(i) quan- 
tifies the degree of synergy present and is /3d-r(i) = 
in the absence of synergistic effects when the model re- 
duces to the simple SIR process. The expression for 
/3d-r(i) depends on the type of synergy. For r-synergy 
(Fig. [IJa)), we assume that /3d-i{t) = Piririt) ~ 1), 
where nr{t) is the number of neighbours challenging 
a recipient host. The rate /3 gives an effective mea- 
sure of the strength of synergy which is constructive 
for /3 > and interfering if /3 < 0. Such a form of 
Pd-r(t) ensures that Ad-r = a for non-synergistic trans- 
missions with rir = 1. For d-synergy (Fig. Hl^b)), we 
assume that /3d-r(i) = f3nd(t), where nd{t) is the num- 
ber of infected neighbours connected to a donor at time 
t so that Ad-r = C( for an isolated infectious host with 
rid = 0. 

Fig. [5] shows the phase diagram for epidemics start- 
ing from a single infected host placed in the centre 
of the lattice. The spread of infection has been nu- 
merically simulated by a continuous-time algorithm 
which is an extension of the n-fold way algorithnt^S. 
The threshold for invasion defines a line of critical 
points etc (/3) separating the non-invasive regime where 
the probability of invasion is Pinv = (in an infi- 
nite system) from the invasive regime characterised by 
-Pinv > 0^ . The phase boundaries shown in Fig. [2] 
correspond to the values of q;c(/3) for i — > oo. The 
finite-size effects were accounted for and eliminated by 
means of finite-size scaling (sec Appendices [XI and iB)) . 

For both types of synergy, ac is a non-increasing 
function of /3, as expected from the monotonic depen- 
dence of Ad-r on l3 [Eq. ([1])]. In the absence of synergy, 
the invasion thresholds for both types of synergy co- 
incide with ac(/3 = 0) = r"^ In 2 (cf. Fig. [2] and Ap- 
pendix [A|. The larger deviations of ac{(3) from ac(0) 
for a given /3 observed for d-synergy are due to the fact 
that, except for the initially infected host, d-synergy is 
present in every transmission event with (3^0 because 
rid > 1 for at least some time. In contrast, for r-synergy 
to be operative in a transmission event, there must be 
rij. > 2 attacking neighbours which is not necessarily 
the case in every transmission event. 

For r-synergy with large positive values of /S, the 
critical line Q!c(/3) tends towards the limiting value 
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FIG. 2. Phase diagram for synergistic epidemics. The ex- 
act phase boundaries are marked by squares for r-synergy 
and by circles for d-synergy. Epidemics are invasive if the 
parameters (a, /3) are on the right of the phase boundary. 
The arrow indicates the limiting value of Oc — 0.2 reached 
asymptotically by the phase boundary for r-synergy in the 
limit of large /3 (not shown on the scale of the graph). The 
approximate phase boundaries obtained analytically by ne- 
glecting correlations in transmission of infection (cf. Ap- 
pendix [Fj are shown for both r- (solid line) and d-synergy 
(dashed line). 



ac(oo) ~ 0.2 (Fig. [2|). In this situation, Ad-r = a < oo 
for rir = 1 but Ad-r = oo if rir > 1, meaning that 
hosts being simultaneously attacked by more than one 
neighbour are infected immediately. In cases with large 
interference in transmission (i.e. very negative /3), the 
invasion threshold is located at ac = 0.86 ± 0.01 inde- 
pendently of the value of /3. In this regime, Ad-r does 
not depend on f3 and corresponds to the limiting sit- 
uation with Ad-r = a for rir = 1 and Ad-r = for any 
rir > 1. 

For d-synergy and values of /3 > 0.8, invasion is pos- 
sible for any positive a. The condition a > is nec- 
essary in order for the epidemic to start from a single 
inoculated site. Once the pathogen is transmitted to 
one of the neighbours of the initially inoculated host, 
rid = 1 for the newly infected host. The combination 
of rid = 1 and synergy is sufficient to make invasion 
possible, irrespective of the value of a. 

Systems with constructive synergy are characterised 
by dense patterns of invasion (Fig. [3ja)). In contrast, 
the invaded region is more sparse for interfering syn- 
ergy (Fig. [3l^b)). This scenario applies to both d- and 
r-synergy. Despite the fact that patterns for interfering 
synergy arc more sparse, it is striking that epidemics 
with interfering synergy can nevertheless be as inva- 
sive as epidemics with constructive synergy in terms 
of the spatial extent of infection. For instance, all the 
patterns shown in Fig. [3] have the same probability of 
invasion, Pjnv = 0.5. In terms of the mean spatial den- 
sity of invasion defined as the relative number of hosts 
infected by the pathogen before invasion occurs with a 
given probability, we conclude that the larger the inter- 
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ferencc in transmission of the pathogen, the less dam- 
aging the invasion is (Appendix [C]). This result can be 
qualitatively understood as follows. Interfering synergy 
favours transmission of infection to hosts with few in- 
fected neighbours and disfavours transmission to hosts 
with several previously infected neighbours. Therefore, 
infection has a tendency to evolve towards poorly in- 
fected regions rather than infecting as many hosts as 
possible. These mechanisms are qualitatively similar 
for both types of synergy but the density of infection 
is larger for d-synergy than for r-synergy for any P ^ Q 
(see Appendix [C] for more detail) . 

For constructive synergy, the patterns of invasion are 
not very much influenced by the particular type of syn- 
ergy. For interfering synergy, the degree of branching of 
the paths followed by the pathogen is clearly smaller for 
d-synergy than for r-synergy (cf. panels in Fig. [Sfb)). 
Branching is always possible for r-synergy. In contrast, 
the patterns of invasion for d-synergy display a branch- 
ing transition for a value of a = Q!b(/3) = ^P'- branch- 
ing occurs for a > ab(/3) but it is absent for a < Q;b(/3) 
(Fig. |4|a)). In the later case, the trajectories of inva- 
sion are of the type followed by a growing self-avoiding 
walk (SAW) ^^'-^^ with an example shown in Fig. SJ^b). 
As expected for growing SAWs, the pathogen can dis- 
play self-trapping, meaning that the epidemic stops if 
the infection reaches a host surrounded by hosts that 
have already been infected (Fig. IDJb)). For values of 
a > ahiP) (Fig. EJb), lower panel), Ad-r = a -f- /3 is 
positive for = 1 and branching is possible. 

As a measure of the temporal ejficiency for invasion, 
we consider the time, ^inv, it takes for the pathogen 
to invade the system. Numerical simulations (see Ap- 
pendix [D|) show that for any given value of Pinv, the 
time tinv decreases with increasing /? (for both types of 
synergy). Therefore, systems with interfering synergy 
are less time-efficient than those with constructive syn- 
ergy. The largest deviations of ijnv from its value with- 
out synergy arc obtained for epidemics with d-synergy 
that operates in every transmission event. 

The analysis presented above demonstrates that syn- 
ergy in transmission of infection has significant and 
sometimes paradoxical and unexpected effects on epi- 
demics. Despite the simple assumptions of the model 
(such as, e.g., short-range synergy and linear depen- 
dence of the rate of infection on the number of in- 
fected nearest neighbours), it reproduces explorative 
and exploitative foraging strategies that are typically 
observed in bacterial, fungal, and tumour growth^r^i. 
The explorative/exploitative behaviour in our model is 
linked to interfering/constructive synergy. The forag- 
ing strategy adopted by fungi, bacteria or ants is known 
to be explorative/exploitative when resources are lim- 
ited/abimdant. 

Changes in the foraging strategy are ultimately due 
to the spatial correlations in transmission rates emerg- 
ing as a consequence of synergy. Spatial correlations 
in Ad-r appear because the neighbourhood of suffi- 
ciently close pairs of hosts have common nodes and 
thus the rates for each pair are not mutually indepen- 



a Constructive (/3 > 0) \) Interfering (/3 < 0) 




FIG. 3. Patterns of infection. Illustration of the effect 
of (a) constructive and (b) interfering synergy on the pat- 
terns of infection spreading from the central host (marked 
by solid circle) in systems of linear size L = 31. The synergy 
rates are /3 = 5 and P = —5 for patterns with constructive 
and interfering synergy, respectively. All snapshots corre- 
spond to the final state of the epidemic with only R and S 
hosts remaining. Solid lines indicate those edges between 
d-r pairs that have transmitted the pathogen at some time 
during the course of the epidemic. The value of a has been 
chosen in each case such that the probability of invasion 
IS -finv = 0.5 for all snapshots: (a) a = 0.47 for r-synergy 
and Q = 0.18 for d-synergy. (b) a = 0.81 for r-synergy and 
a — 5.22 for d-synergy. 



dent. Synergistic epidemics can then be regarded as 
a correlated dynamical percolation analogous to the 
well-known mapping of the non-synergistic SIR pro- 
cess to uncorrelated dynamical percolation (see— i^ii^>^ 
and details on the mapping to correlated dynami- 
cal percolation in Appendix [E|. In most situations, 
spatial correlations in transmission are short-ranged 
and the critical behaviour of epidemics at etc belongs 
to the dynamical uncorrelated bond-percolation uni- 
versality clas a^'^^ . However, for large interfering d- 
synergies with a < ab(/3) (region under continuous 
line in Fig. IDJa)) correlations become effectively long- 
ranged and epidemics are growing SAWs whose critical 
properties belong to the universality class of the stan- 
dard SAW^^. Although the SAW behaviour affects 
the local properties of epidemics with a > ab(/3), the 
large-scale behaviour at ac is not affected (i.e. the 
critical exponents at invasion are the same as those for 
uncorrelated bond percolation, as shown in detail in 
Appendix iBjl . This is a consequence of the fact that 
the probability that a growing SAW invades a large 
system is known to be zero due to the self-trapping 
phcnomcno n^^d^ . This implies that Uc > for any 
value of /3 (cf. Fig.gKa)). 
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FIG. 4. Branching and invasion transitions for d-synergy. 
(a) Diagram showing the branching transition hne ab(/3) = 
— (3 together with the phase boundary for invasion/non- 
invasion transition, ac(/3)- Branching is forbidden in the 
region under the continuous line where a < Qc(/3)- In- 
vasive epidemics can occur in the region above the dashed 
line corresponding to a > ac(/3) where branching is present. 
The intermediate region between the continuous and dashed 
lines with ab < a < ttc corresponds to epidemics that dis- 
play branching but are not invasive, (b) Example of path 
for infection with d-synergy in the regime without branch- 
ing {a < ah). Infection starts from the host marked by a 
solid circle and evolves along the path indicated by arrows. 
Arrow numbers define schematically a sequence of infection 
events. For the first infected host (solid circle), rid = and 
infection can be transmitted to any of its neighbours at a 
rate Ad-r = a. However, rid = 1 as soon as the pathogen 
is transmitted to one of the neighbours (arrow 1) and thus 
Ad-r ~ because a < —j3. At this moment, transmission 
of infection is arrested until the initially infected host re- 
covers. After this recovery, rid = for the newly infected 
host and infection can be transmitted to one of its nearest 
neighbours at rate a. Iteration of this process over time 
gives growing SAWs. The state after the event marked by 
arrow 10 illustrates the phenomenon of self-trapping. 



The proposed model becomes analytically tractable 
if spatial correlations in transmission are assumed to be 
negligible. The analytic solutions provide a good qual- 
itative description of the main features of the phase 
diagram for both r- and d-synergy. Correlations in the 
exact model prevent the approximate description from 
being quantitative (cf . Appendix |F] for a complete de- 
scription). 

In summary, the presented work shows that synergis- 
tic effects at the individual level play an important role 
in invasion at the population level. The analysis has 
been restricted to the spread of epidemics in 2D regu- 
lar networks relevant for, e.g., populations of plants in 
a field. The extension of our work to synergistic effects 
for spreading processes in higher-dimensional lattices 
or more complex networks is not only conceptually ap- 
pealing but also important for problems in multiple dis- 
ciplines^. Although such an extension is technically 
straightforward, understanding the interplay between 
synergistic effects and topological properties for differ- 
ent types of networks is a challenging task for future 
work. 
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APPENDICES 
Appendix A: Probability of invasion 

In this section, we discuss the dependence of the 
probability of invasion, Pjnv, on the rates a and j3. The 
probability of invasion is defined as the relative num- 
ber of invasive events out of many (in our simulations, 
> 5000) stochastic realizations of epidemics. Although 
we will deal with systems of finite size, the trends of 
Pinv (see Fig.[S]) are in qualitative agreement with those 
expected from the phase diagram for epidemics in infi- 
nite systems shown in the main text (Fig. 2). For both 
r- and d-synergy, the invasion curves (Pinv '<^s a) ap- 
proach zero for small transition rates a thus describing 
the non-invasive regime for epidemics. For large values 
of a, the invasion probability is finite so that it de- 
scribes the invasive regime. The invasion curves start 
to deviate from zero at progressively lower values of a 
as /? increases. This illustrates the non- increasing char- 
acter of q;c(/3) for increasing (3 discussed in the main 
text. The tendencies of Pi„v in extreme situations re- 
quire further explanations that are different for r- and 
d-synergy. 

a. r-synergy Fig. [5ja) shows the probability of 
invasion for epidemics with r-synergy in systems of lin- 
ear size L = 31 as a function of a (invasion curves) for 
several values of /3. For large positive values of /3, the 
invasion curves tend to a single (master) curve which 
can be described by a limiting function of a that does 
not depend significantly on /3. This is illustrated in 
Fig- Eta) by two curves for /3 = 50 and (3 — 100 coincid- 
ing within numerical error. Consequently, on the phase 
diagram, the critical line ac(/3) tends towards the lim- 
iting value ac(oo) — 0.2 (shown by the arrow in Fig. 2, 
main text). In cases with large interference in trans- 
mission (corresponding to very negative /3) , the curves 
for Pi,iv(Q!) again collapse on a limiting curve shown in 
Fig. El^a) for (3 = —1.5 and j3 = —2.0. The correspond- 
ing invasion threshold in the thermodynamic limit is 
located at ac = 0.86 ± 0.01 independently of the value 
of /3. 

h. d-synergy For d-synergy and values of /3 > 
0.8, invasion is possible for any positive value of a. In 
this regime, the invasion curves collapse on a master 
curve that does not depend on the value of /? (cf. curves 
in Fig.[51[b) for /3 = 5 and /3 = 50). For negative values 
of /?, the invasion curves do not approach a limiting 
function as /3 decreases (Fig. [S^b)). As a result, Uc 
increases monotonically with decreasing j3 (see Fig. 2, 
main text). 
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FIG. 5. Probability of invasion Pinv in a system of linear size L = 31 as a function of the elementary rate a for (a) 
r-synergy and (b) d-synergy. Different curves correspond to different values of the synergy rate, /3, as indicated in the 
legend for each case. 



Appendix B: Critical behaviour and finite-size 
effects 

In this section, we give details about the methods 
used for obtaining the phase diagram for invasion in the 
limit L — > c» (sec Fig. 2 in the main text). In addition, 
we give numerical support for the statements made in 
the main text about the critical behaviour displayed 
by epidemics with d-synergy in the different regions 
presented in Fig. 4(a) of the main text. 

The data points for the phase boundary, adP), 
in the phase diagram for invasion were obtained by 
analysing both the probability of invasion, Pinv, and 
the average relative number, A^i, of epidemics spanning 
the system in one and only one direction (ID-spanning 
epidemics). The analysis of A^i presented here is sim- 
ilar to that performed for avalanches in spin models 
with quenched disorde r^^i^^ . 

Below we show that, for epidemics with constructive 
or weakly interfering synergy, the method based on the 
analysis of ID-spanning clusters gives more accurate 
estimates for adP) than those obtained by using Pinv 
In particular, our estimates for the invasion threshold 
for non-synergistic epidemics based on A'^i are in ex- 
cellent agreement with the value ac(0) = r~^ln2 ex- 
pected from the mapping to bond-percolation valid for 
(3 = O'^ . In contrast, for epidemics with large interfer- 
ing d-synergy, the estimation of ac {l3) based on Pinv is 
more accurate than that achieved by using A^i. 



1. Epidemics with constructive or wealtly 
interfering synergy 

For epidemics with constructive or weakly interfering 
synergy, the quantity A^i has a peak when plotted as a 
function of a for a fixed value of /?. This is the case for 
both r- and d-synergy. Fig. [S] shows a representative 
example corresponding to constructive d-synergy with 
/3 = 0.3. The position of the peak for A^i is related to 
the critical value of the parameter a separating invasive 



and non-invasive regimes. Indeed, in the non-invasive 
regime, TVi ~ because the clusters of removed hosts 
are of small size and thus they have a negligible prob- 
ability to touch the boundaries of a finite system. On 
the other hand, iVi is expected to be close to zero in the 
invasive regime as well because the pathogen spreads 
in two directions rather than along one only. As a 
result, A'^i can be finite at the invasion threshold or, 
due to finite-size effects, in a certain region around the 
threshold. The critical value q;c(/3) in the limit L — > oo 
is obtained using the scaling properties of A^i(a, L) in 
the vicinity of the invasion threshold. Due to the ab- 
sence of characteristic length scales at criticality, the 
dependence of Ni on a and L is expected to obey the 
following scaling law^i*^: 

iVi(a,i) =i^iVi((a-ae)ii/") , (Bl) 

where 9 and v are critical exponents and iVi is a scaling 
function which depends on a and L through the prod- 
uct {a — ac)L^^'^ only. The values of 9, v, and ac can 
be determined by scaling collapse for Ni{a,L). Tech- 
nically, this can be achieved by plotting the quantity 
L~^Ni vs (a — ac)L^^^ and requiring that the scal- 
ing hypothesis (jBl[) is satisfied (i.e. curves for differ- 
ent L collapse on a single master curve corresponding 
to A^i). Fig. [61(b) shows an example of such scaling 
collapse for d-synergy with /3 = 0.3. In this case, 
the collapse gives the values ac = 0.445 ± 0.005 and 
9 = —0.10 ±0.05. Bearing in mind the mapping of SIR 
epidemics to dynamical uncorrclated bond-percolation 
holding for /3 = 0^, the value of v has been fixed to 
V = A/i corresponding to uncorrclated percolatiori^. 
The good quality of the collapse suggests that, despite 
the existence of correlations in transmission of infec- 
tion for /3 7^ 0, the critical behaviour at q;c(/3) falls into 
the universality class for the uncorrclated percolation. 
This is the expected behaviour when correlations are 
short-ranged^. 

The above statement is also supported by the be- 
haviour of Pinv in the vicinity of the invasion threshold 



6 




0.15 



CD 




0.05 



(a-a„) L 



FIG. 6. (a) Dependence of the average relative number of ID-spanning epidemics, TVi, on the elementary rate a for d- 
synergy with /3 = 0.3. Different curves correspond to systems of different size L, as marked in the legend, (b) Scaling 
collapse of A'^i according to the scaling law (|Bip for the curves shown in (a). The critical exponents have been set to the 
values v — 4/3 and /3 = 5/36, corresponding to values found for uncorrelated percolation. The scaling collapse gives the 
invasion threshold Qc — 0.445 ± 0.005. 
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FIG. 7. (a) Dependence of the probability of on the elementary rate a for d-synergy with /3 = 0.3. Different 

curves correspond to systems of different size L, as marked in the legend, (b) Scaling collapse of Pinv according to the 
scaling law given by Eq. (|B2|) for the curves shown in (a). The critical exponents have been set to the values v = 4/3 and 
P = 5/36. The invasion threshold has been set to the value ac = 0.445 obtained from the collapse of A''i (cf. Fig. [6]). 



which satisfies the following sealing hypothesis, 

Pinv = L-^/'^Pin.iia - «c)i'/"), (B2) 

with the exponents = 4/3 and (3 — 5/36 eorre- 
sponding to uncorrelated percolation. As an exam- 
ple, Fig. [TJb) shows the scaling collapse for the curves 
shown in Fig. [Tlja) for several system sizes. The col- 
lapse in Fig. [Tl^b) has been obtained by setting ac to 
the value ac = 0.445 obtained from the collapse of A^i. 
The quality of this collapse is remarkably good despite 
the fact that no free parameters have been used (i.e. 
the value of i^, /3, and ac has been considered as being 
fixed). 

If instead of fixing a^, we try to estimate its value 
from the collapse of Pinv, the sigmoidal shape of Pinv 
prevents the estimate from being as accurate as the 
one obtained from the collapse of Ni{a,L) which are 
peak-shaped. 



2. Epidemics with large interfering d-synergy 

In the main text, we have shown that epidemics with 
interfering d-synergy behave as growing self-avoiding 
walks (SAWs) for a < a^ = —(3 (cf. Fig. 4(a), main 
text). Fig. [Ha) shows Ni vs a for epidemics with 
/3 = — 10 spreading in systems of different size. As can 
be seen, growing SAWs give a constant contribution to 
iVi for a < 10 = ab(/3). The probability that a grow- 
ing SAW spans the system in one dimension decreases 
with the system size, L, and thus the contribution of 
such objects to Ni tends to zero as i — > oo. In our sim- 
ulations, the systems are finite and contribution of the 
ID SAWs to Ni is not negligible. As a consequence, 
the scaling collapses of A^i based on the hypothesis by 
Eq. (jBl[) are of low quality and not very useful for the 
estimation of Q!c(/?) for some negative values of /?. 

In this situation, the estimates of the critical val- 
ues of the rate a are more conveniently obtained by 
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FIG. 8. (a) Dependence of the average relative number of ID-spanning epidemics, A''i, on the elementary rate a for d- 
synergy with P = —10.0. Different curves correspond to systems of different size L, as marked in the legend, (b) Scaling 
collapse of A^^i according to the scaling law given by Eq. (|Bl| for the curves shown in (a). The critical exponents have been 
set to the values — 4/3 and /3 = 5/36. The invasion threshold has been set to the value cic ~ 10.15 obtained from the 
collapse of Pinv (cf. Fig. ^ . 
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FIG. 9. (a) Dependence of the probability of on the elementary rate a for d-synergy with /? = —10. Different 

curves correspond to systems of different size L, as marked in the legend, (b) Scaling collapse of Pinv according to the 
scaling law given by Eq. (|B2ll for the curves shown in (a). The critical exponents have been set to the values v = 4/3 and 
/3 — 5/36, corresponding to uncorrelated percolation. The scahng collapse gives the invasion threshold Qc = 10.15 ± 0.03. 



analysing the probability of invasion, Pinv. This quan- 
tity is less afFected by the growing SAW epidemics be- 
cause invasion requires reaching all the four edges of the 
system. Fig.[9l[a) shows the dependence of Pi„v on a for 
/3 = —10. Fig. [9jb) demonstrates the scaling collapse 
of the curves displayed in panel (a) with the exponents 
= 4/3 and (3 = 5/36 corresponding to uncorrelated 
percolation and etc = 10.15 ± 0.03. The quality of the 
collapse suggests that, despite the proximity of Uc to 
ttb, the behaviour at large scales corresponds to that of 
dynamical uncorrelated percolation. Use of the same 
value of the critical rate, ac = 10.15, for Ni{a,L) re- 
sults in a collapse of much poor quality (see Fig. [HI)- As 
expected, the quality of the collapse is reasonable only 
for relatively large values of (a — ac)L^^'^ for which the 
influence of growing SAWs is negligible. 



Appendix C: Mean spatial density of invasion 

In this section, we give further quantitative support 
to the results discussed in the main text concerning 
the effect of synergy on the spatial density of invasion. 
The mean spatial density of invasion is defined as the 
relative number of hosts, tir, that are in the removed 
state (R) by the end of an invasive epidemic (i.e. the 
relative number of hosts that have been infected dur- 
ing the course of the epidemic and are removed by the 
end) . Due to stochasticity in the transmission of infec- 
tion, different realisations of invasive epidemics charac- 
terised by the same parameters a and /3 have a different 
random value for np . The probability density function 
for the density of infection, p{n^), has a single peak for 
any values of the parameters a and P (see the inset in 
Fig. [T0|) . Therefore, riR has a well-defined scale that, 
due to the small degree of asymmetry of p(n^), can be 
properly represented by the mean (rig) with dispersion 



8 



s 




-6 -4 -2 2 4 6 8 10 -6 -4 -2 2 4 6 8 10 



FIG. 10. Mean density of invasion, tir, as a function of 
P for epidemics with Pinv = 0.5 in systems of linear size 
L = 31. Symbols (squares for r-synergy and circles for 
d-synergy) correspond to the mean of wr over stochastic 
realisations and error bars show the standard deviation. 
The inset displays representative examples of the probabil- 
ity density function for the mean density of invasion, p{nj\), 
for P = -2. 



FIG. 11. Time to invasion, tinv, as a function of /3 for 
epidemics with Pinv = 0.5 in systems of linear size L = 31. 
Symbols (squares for r-synergy and circles for d-synergy) 
correspond to the mean of finv and error bars show the 
standard deviation. The inset displays representative ex- 
amples of the probability density function for the time to 
invasion, pitinv), for j3 = —2. 



Appendix D: Temporal efficiency 



given by the standard deviation. Fig. [TO] shows the 
dependence of the spatial efheiency on (3 for invasive 
epidemics with value of a chosen in each case so that 
Pinv — 0.5. As stated in the main text, the mean den- 
sity of invasion exhibits a global tendency to increase 
with increasing (3 both for r- and d-synergy. 

As expected, is identical for the two types of syn- 
ergy if /3 = 0. For any non-zero value of /3, the density 
of infection for d-synergy, riR^dj is larger than that for 
r-synergy, nR^i-. In cases with constructive synergy, it 
is likely that riR.d > 'T'R,!- because, as argued in the 
main text, synergistic effects are more prominent for 
d-synergy (they operate in every transmission event). 
A plausible explanation for the origin of the inequality 
'^R.d > J^R,!- for /? < can be given by recalling that 
branching is more frequent in paths of infection for r- 
synergy than for d-synergy. Due to the higher degree of 
branching for r-synergy, it is more probable that hosts 
become effectively isolated from infection for this type 
of synergy if they are simultaneously challenged by sev- 
eral neighbours which interfere and do not transmit 
infection. As an extreme case, consider a situation in 
which a host is simultaneously challenged by its all four 
neighbours. It is clear that the challenged host will be- 
come inaccessible forever if infection is not transmitted 
by any of the four challenging neighbours. In contrast, 
the lower degree of branching for d-synergy makes the 
existence of effectively isolated hosts less likely. As a 
consequence, more hosts can be infected for d-synergy 
during the course of epidemics. 



This section complements the part of the main text 
devoted to the time to invasion, Unv- Due to stochas- 
ticity in the transmission of infection, t-mv is a random 
variable described by a probability density function, 
^(iinv), which has a single peak for any values of the 
parameters a and j3 (see the inset in Fig. Ilip . We can 
then proceed in an analogous manner as we have done 
above for the mean density of infection and describe 
iinv by its mean (tinv) and dispersion given by the stan- 
dard deviation. Fig. [11] shows that for both types of 
synergy (^inv) decreases with increasing /?. The dis- 
persion of p(tinv) also decreases with increasing /3, as 
indicated by the error bars in Fig. [TTl 

By definition, the time to invasion is identical for the 
two types of synergy if /3 = 0. For /3 ^ 0, the largest de- 
viations of <inv from its value without synergy (/? = 0) 
(cf . circles and squares in Fig. [TT|) are for epidemics ex- 
hibiting d-synergy. This is mostly due to the fact that 
d-synergy operates in every transmission event. In av- 
erage, this makes transmission events slower/quicker 
for interfering/constructive d-synergy. In addition to 
this factor, it is likely that the higher degree of branch- 
ing in the foraging strategy for infection with interfer- 
ing r-synergy also contributes to making <inv smaller 
for r-synergy with /? < 0. 



Appendix E: Invasion as a correlated dynamical 
bond-percolation problem 

The aim of this section is twofold. First, we give de- 
tailed definitions for the transmissibility in synergistic 
epidemics and related quantities such as its probabil- 
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ity density function (p.d.f.) and mean value. Second, 
the effect of correlations on the probability of invasion 
summarised in the main text is analysed here in more 
detail. 



1. Transmissibility 

In the presence of synergistic effects, the trans- 
mission of infection can be described as a non- 
homogeneous Poisson process with the time-dependent 
infection rate \d-i{t) defined by Eq. (1) in the main 
text. The probability that the infection has not been 
transmitted in a d-r (donor-recipient) pair by time t 
defines the survival probability Sd-rit). For a non- 
homogeneous Poisson process, Sd-i{t) satisfies the fol- 
lowing differential equatio n^^'^^ : 



d5, 



d-r 



dt 



= -Ad-r(t)5d-r(0 



The solution of this equation with initial condition 
S{t = 0) = 1 (such a condition ensures that the 
pathogen is not transmitted instantaneously when the 
donor is infected) is 



Sd-rit) = exp 



- / Ad-r(i) dt 



(El) 



The transmissibility T^.r is defined as the probability 
that the pathogen is transmitted from the donor to 
the recipient over the infectious period of the donor, r. 
Therefore, it can be expressed in terms of the survival 
probability as follows: 



l-^d-r(T) . 



(E2) 



Substitution of the expression for Sd-r given by 
Eq. (|Eip into Eq. (jE2|l results in the following expres- 
sion for Td-i '. 



Td 



1 — exp 



Ad-r(t) dt 



(E3) 



In the absence of synergistic effects, the rate 
Ad-r(t) = « remains constant over the infectious pe- 
riod T and is the same for all d-r pairs. In this case, 
the transmissibility reduces to the homogeneous value, 
Td-r = 1— e""^", that plays a central role in the mapping 
of SIR epidemics to the well-known uncorrelated dy- 
namical percolation^ii^. In this mapping, Td-r is iden- 
tified with the bond probability and Pinv from an ini- 
tially inoculated site is identified with the probability 
Poo that such site belongs to the infinite cluster of con- 
nected sites in the dynamical percolation problem^. 
For SIR epidemics, Pinv is fnlly parameterised by Td-r. 
This is analogous to the fact that Poo is fully parame- 
terised by the bond probability in dynamical percola- 
tion. 

When synergy is present, the transmissibility Td-r is 
a functional that depends on the time evolution of the 



rate \d-i{t) which, in turn, depends on the infection 
history of the neighbouring hosts to the d-r pair. Each 
d-r pair involved in an epidemic is in general charac- 
terised by a different dependence of the rate Xd-r{t) 
on time. As a consequence, the field of transmissibil- 
ities is spatially heterogeneous. In addition, the de- 
pendence of Td-r on the neighbourhood of the d-r pair 
introduces non-trivial correlations in transmission and 
thus in transmissibilities. 

In order to study the role of synergy-induced corre- 
lated heterogeneity at the host level on Pinv we proceed 
in a way inspired from previous works dealing with het- 
erogeneous SIR epidemicsii^'^ . The simplest situ- 
ation with heterogeneity in transmission corresponds 
to epidemics where Ti-r are independent random vari- 
ables for all d-r pairs. In this case, Pinv only depends 
on the mean transmissibility, (T]^"^. In more com- 
plicated situations, the transmissibilities for different 
d-r pairs are not independent and Pinv is a function of 
the whole set of transmissibilities, {Td-r}, that cannot 
be completely parametrised by (r)ii^^ An exact 
mapping of such SIR epidemics to uncorrelated per- 
colation is not possible in general. However, the use 
of (T) can still be useful to analyse the consequences 
that heterogeneity in local transmission has on Pinv 
For instance, a considerable progress has been made 
in understanding the role of correlations in epidemics 
where heterogeneity in transmission is associated with 
heterogeneity in recovery times of infected hosts. In 
this case, the following important result has been rig- 
orously derived^!"— : for a given value of (T), the re- 
silience to invasion increases with increasing degree of 
heterogeneity (more precisely, P^^^{{T)) < PtnT(iT)), 
where and Pinv™ a^re the probabilities of invasion 
for heterogeneously and homogeneously distributed re- 
moval times, respectively). Here, we show that the 
dependence of Pinv on (T) for synergistic epidemics is 
more complicated (see Fig. [T2|) . In spite of that, in 
subsections IE 21 and IE 31 we show that analysing the de- 
pendence of Pinv on the degree of synergy for given (T) 
is still informative. 

For synergistic epidemics we define the mean trans- 
missibility in terms of two averages: spatial average in 
each realisation of epidemics and stochastic over dif- 
ferent epidemic realisations. The spatial average for 
transmissibility, T^, is calculated for a particular r-th 
realisation of the epidemic in the following manner. 



1 



Nd- 



d-r 



where the sum extends over the number iVd-r of d-r 
pairs in the final state of the epidemic, i.e. over all d-r 
pairs challenged by the infection. The value of Td-r in 
the above equation for each d-r pair is calculated using 
Eq. (|E3|) with the transmission rate Xd-r{t) measured 
numerically for the r-th realisation of epidemic. The 
direct numerical evaluation of the transmissibility for 
synergistic epidemics as a frequency of successful trans- 
mission of infection between donor and recipient in the 
d-r pair would require reproduction of time-dependent 



10 



transmission rates giving exactly the same integral over 
time, — Xd-r{t)'it (i.e. all possible rates Xd-i {t) giv- 
ing the same value of Td-r in Eq. (jE3[l ). This imposes 
a very restrictive condition on the time of infection of 
the nodes in the neighbourhood of the d-r pair and 
thus evaluation of T^.^ as a frequency can be hardly a 
feasible computational task. The mean transmissibil- 
ity, (T), is obtained by stochastic averaging of Tr for 
R different realisations of the epidemic: 

r—1 

where the value of Tr is averaged over R differ- 
ent stochastic realisations of the epidemic. Note 
that averaging over stochastic realisations is necessary 
to account for the fact that different realisations of 
epidemics lead to different spatial configurations for 

{7d-r}. 

Bearing in mind that the population of hosts is ho- 
mogeneous (i.e. a and /3 do not depend on the d-r pair 
location) the values of transmissibilities for any d-r pair 
given by Eq. (|E3p are independent random variables 
characterized by a p.d.f., p{Td-r)- Once p(Td_r) is avail- 
able, e.g. numerically, the mean transmissibility can 
equivalently be calculated as 

(T) = / Tp{T) dT . (E4) 
Jo 

The analysis of p(rd-i ) is important for understanding 
the dependence of Pinv on the degree of synergy and 
(T), as we show below for r- and d-synergy. 

2. r-synergy 

Fig. [TH shows the dependence of Pjnv on (T) for 
several values of p. For any fixed value of (T), the 
probability of invasion with interfering synergy (see the 
curves corresponding to /3 < and marked by open and 
solid circles and open squares in Fig. [T2|) is systemat- 
ically greater than for non-synergistic epidemics (see 
the curve for /3 = marked by solid circles in Fig. [T2|) . 
i.e. the populations exhibiting interfering r-synergy 
are more vulnerable to invasion than those without 
synergy. For constructive r-synergy, it is possible to 
distinguish between two different regimes. The first 
regime corresponds to epidemics with moderate syn- 
ergy (curves for (3 = 0.8 and 2 in Fig. [12]) that are less 
invasive than non-synergistic epidemics for all values of 
(T), i.e. the curves marked by the solid diamonds and 
stars are always below the curve marked by the solid 
squares. The second regime corresponds to larger val- 
ues of (3 (e.g. /3 = 50 and 100 in Fig. [T^ . In this case, 
there is a range of (T) where the curves marked by 
open triangles and crosses are above the curve marked 
by the solid squares. In this interval for (T), the epi- 
demics with constructive r-synergy are more invasive 
than those without synergy. 




<T> 

FIG. 12. Probability of invasion as a function of the mean 
transmissibility for r-synergy in systems of linear size L — 
31. Different curves correspond to different values of /? as 
marked in the figure. 



These features can be qualitatively illustrated by 
analysing the evolution of the distribution of trans- 
missibilities, p(rd-r), with the strength of synergy. 
Without synergy, p(T(j_i.) has a (5-functional shape, 
p{Td-r) = S{Td., - Ti), where Ti = 1 - e""^ is the 
non-synergistic transmissibility. The epidemic is inva- 
sive (non-invasive) if Ti > = 1/2 (Ti < T^)^. Once 
r-synergy is introduced, the shape of p{Td-i ) changes. 
The (5-functional peak is still present and it describes 
recipients with a single infected neighbour, nr = 1. 
In addition, new contributions to /9(Td-r) appear for 
Td_i. < Ti (Td_r > ^i) in case of negative (positive) 
values of /3 (see the lower panel in Fig. [13] where the 
(5-functional peak at Td-r = 7i is marked by o). Such 
contributions come from recipients surrounded by more 
than one infected neighbour, i.e. tIj. > 1. The num- 
ber of infected neighbours varies in discrete manner 
and this brings a non-smooth functional dependence 
to /9(Td-r) consisting of cusps associated with the dis- 
crete changes in ny(t) and smooth components between 
cusps originated from the continuity of time. In other 
words, nj-(t) is a piece- wise function of time the integral 
of which produces a continuous set of transmissibilities 
according to Eq. (|E3|1 . 

For interfering synergy, the average transmissibil- 
ity is (T) < Ti, due to the contribution to p(Td-r) 
of Td_r < Ti described above. If the heterogeneous 
transmissions associated with synergy were uncorre- 
lated, the probability of invasion for synergistic epi- 
demics with the elementary rate a would be smaller 
than that corresponding to epidemics with the same 
value of a but without synergy. This is a consequence 
of the fact that the probability of invasion for systems 
with heterogeneous but uncorrelated transmissions de- 
pends on (T) only^i^i^. However, this is not the case 
for epidemics with interfering synergy which are more 
invasive than epidemics without synergy. This is due to 
correlations between synergistic transmissibilities. In- 
deed, the d-r pairs that passed the infection in an epi- 
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FIG. 13. The upper panels show Pi,w vs. (T) for four values of /3 for r-synergy. The lower panels show the p.d.f. 
p(rd-r) corresponding to the points in the invasion curves indicated in the upper panels with an arrow. The mean of the 
p.d.f. p(rd-r) gives the value of (T) in the upper panel (Eq. (|E4[) \ The non-synergistic transmissibility Ti = 1 — e"""^ 
corresponding to the recipients with rir = 1 is marked by o in all the panels. The critical transmissibility for non-synergistic 
epidemics is Tc — 1/2. The intervals of transmissibility with T < Tc and T > Tc are indicated by the blue and orange 
shaded regions, respectively. 



demic (mostly those with non-synergistic transmissibil- 
ity Ti) arc arranged in a spatially correlated finger- like 
manner (as shown in Fig. 3(b) in the main text) that 
makes invasion possible. Therefore, the synergistic epi- 
demics can be mapped onto the correlated dynamical 
bond-percolation problem in which the bond probabil- 
ities are associated with transmissibilities. 

For constructive synergy with moderate values of f3 
(corresponding to the first regime mentioned above), 
most of the d-r pairs in the system have transmissibility 
Ti < Tc = 1/2 (cf. the position of the peak marked by o 
with the rest of the distribution in the panel for /3 = 0.2 
in Fig.flS]). Most of the synergistic d-r pairs have Td-r > 
Tc as can be seen from the comparison of the areas 
under the curve for p.d.f. for Td-r < Tc (excluding non- 
synergistic transmissibilities under the peak marked by 
o) and Td-r > Tc- However, the abundance and value 
of the transmissibilities for such synergistic pairs does 
not seem to be high enough as to allow for invasion 
unless (T) is clearly larger than Tc- 

For larger values of (3, invasion is possible for (T) < 
Tc because, as shown in Fig. [13] for (3 = 5, the relative 
number of synergistic pairs is large enough and they 
have T ~ 1 so that the transmission of infection is very 



likely. Moreover, these pairs are placed in a spatially 
correlated manner which also favours the invasion for 
(T) < Tc. 



3. d-synergy 

The p.d.f. p(Td-r) for epidemics with d-synergy is 
shown in Fig. [M] for several values of (3. At first sight, 
the effects of d-synergy on p{Td-i) are qualitatively sim- 
ilar to those of r-synergy. However, for d-synergy, the 
(5-functional peak for the non-synergistic transmissibil- 
ity To marked by o corresponds to realisations of the 
epidemic when the pathogen has not been transmit- 
ted from the initially infected host to any of its neigh- 
bours (i.e. the epidemic has not started spreading). 
Only in this case nd = over the whole infectious pe- 
riod of the initially inoculated host and Td-r = Tq. In 
contrast, Td-r 7^ Tq for every d-r pair if the epidemic 
starts, meaning that there is no contribution to the 
5-functional peak from these epidemics. 

The scenario for epidemics with interfering d-synergy 
is similar to that for r-synergy, meaning that invasion 
can occur for values of (T) < Tc (cf. the curves for 
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FIG. 14. The upper panels show Pinv vs. (T) for four values of /3 for d-synergy. The lower panels show the p.d.f. 
p(rd-r) corresponding to the points in the invasion curves indicated in the upper panels with an arrow. The mean of the 
p.d.f. p(Td-r) gives the value of (T) in the upper panel (Eq. (|E4[) \ The non-synergistic transmissibility To = 1 — e"""^ 
corresponding to recipients with rid = is marked by o in all the panels. The critical transmissibility for non-synergistic 
epidemics is Tc = 1/2. The intervals of transmissibility with T < Tc and T > Tc are indicated by the blue and orange 
shaded regions, respectively. 



/3 < and ^ = in Fig. [T5|). For any fixed negative 
value of /?, the shift (to the left) of the invasion curve 
from the synergy-free one, Pinv{{T)) with /3 = 0, is 
greater for d-synergy than for r-synergy. Correlations 
in transmission seem to play a very prominent role for 
d-synergy since invasion can occur even for very low 
values of (T). The fact that d-synergy induces larger 
shifts for the Pinv{{T)) curve towards smaller values of 
(T) than r-synergy is due to the greater abundance of 
synergistic connections for d-synergy, as argued in the 
main text. This is clear from the comparison of the 
p.d.f. p{Td-i) for /? < corresponding to the two types 
of synergy shown in Figs. [T3l and [Ml The relative num- 
ber of non-synergistic d-r pairs with non-synergistic 
transmissibility Tq is always smaller for d-synergy. 

For constructive d-synergy, there are two regimes 
that are qualitatively similar to those discussed for con- 
structive r-synergy above. The regime where invasion 
is only possible for (T) > Tc exists for very weak con- 
structive synergy (0 < < 0.1). Figs. [13 and [H] il- 
lustrate the behaviour in this regime for /3 = 0.1. For 
greater values of (3, invasion is possible for (T) < Tc = 
1/2 (see the curves for f3 > 0.4 in Fig. [TS)) . This is 
due to the presence of a large number of synergistic 



transmissions with T^-i > Tc, as illustrated in Fig. [T4| 
for /3 = 0.8. As argued above, the peak at To is due 
to epidemics that do not start spreading. The rest of 
contributions to /9(Td-r) corresponds to those cases in 
which the pathogen is transmitted from the initially 
infected host to at least one of its neighbours. The 
values of the transmissibilities for such epidemics are 
larger and thus Pmv is also larger. In this case, inva- 
sion is possible for any positive value of (T) because 
the epidemic is invasive with high probability provided 
it starts spreading. In other words, many epidemics do 
not start spreading for very low value of a (and thus 
low transmissibility). However, the probability that 
the epidemic starts spreading is larger than zero for 
any a > 0. Once it starts, invasion almost certainly 
occurs. Therefore, invasion is possible for any positive 
value of a, no matter how small the value of a is. 



Appendix F: Phase diagram for a simple model 

The aim of this section is to present a simple model 
for evaluation of the phase diagram in a — /3 parameter 
space for both types of synergy in SIR process. 
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FIG. 15. Probability of invasion as a function of the mean 
transmissibility for d-synergy in systems of linear size L = 
31. Different curves correspond to different values of j3 as 
marked in the figure. 



The phase transition from non-invasive to invasive 
SIR regime in heterogeneous systems with uncorrelated 
transmissibilities occurs when 



(Fl) 



where Tc is a critical topology-dependent bond- 
pcrcolation probability {Tc ~ 1/2 for a square lattice) 
and (T) is the mean transmissibility^^. For a synergis- 
tic SIR process, the transmissibilities are heterogeneous 
due to variable neighbourhood during infection (bond 
creating) process and correlated as it follows from our 
analysis of probability of invasion vs mean transmis- 
sibilities (see Sec. |E|. Such correlations make exact 
analytical treatment for synergistic SIR process hardly 
possible. However, if the correlations between trans- 
missibilities are ignored then a simplified model for 
synergistic SIR process can be introduced and solved 
analytically for the boundaries in the phase diagram. 
The simplified assumptions of the model are the fol- 
lowing: 

(i) there arc no correlations in transmissibilities; 

(ii) the neighbourhood of a r-d pair does not change 
over the infectious period of the donor; 

(iii) the probabilities of various neighbourhoods of a 
d-r do not depend on the rates a and /?; 

Under these assumptions, the analytical expressions 
for the phase boundaries reproducing qualitatively all 
the features found numerically can be derived and anal- 
ysed. 



1. d-synergy 

We start analysis of the simplified model for syner- 
gistic SIR process on a square lattice (for concreteness) 
with the case of d-synergy. Let us consider an infected 
node (other than the initially infected host) which at- 
tempts to transmit infection (create a bond) to one of 
its susceptible neighbours during its infectious period 
r = 1. This process can occur with different probabil- 
ity depending on the number of infected nodes linked 
to the infecting one (donor). Applying assumptions 
(i)-(ii) that transmissibilities are uncorrelated and a 
certain random neighbourhood does not change over 
the infectious period of the donor, we can calculate the 
mean transmissibility in the following way. 



(T) = ^ PnTn 



(F2) 



Here Pn is the probability that the donor is connected 
to n {n = 1,2,3 for square lattice) infected nodes and 
n-depcndent transmissibility, r„, for the d-r pair is 



T = 



1 



-a — 71 [3 







if 
if 
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a - 



> 
n/3 < . 



(F3) 



Different neighbourhoods of the d-r pair occur with 
probabilities p„. Under assumption (iii), these prob- 
abilities do not depend on a and (5 and given by the 
following expressions: 



Pi 

P2 



P3 



3p(l- 


P? 


1-(1- 


-pY 


3P^(1 - 


-p) 


1-(1- 


-pY 








-pY 



(F4) 



where Pn is the probability that the infecting node is 
linked to n nodes given that it is linked at least to one 
out of three possible nodes and p is the bond proba- 
bility with p = Tc at criticality. In Eq. (|F4[) . we used 
the assumption (i) that the bonds were created inde- 
pendently. 

The locus of critical points can be found by solving 
Eq. dni), 



Tc = '^Pn{Tc)Tn , 



(F5) 



n=l 



where (T) is given by Eqs. (IF2])-(|F4|. Eq. ((F5|l can be 

recast as (a > 0), 



- ln(l - Tc) + In {pi{Tc)e~P +p2{Tc)e-^f' + Pz{Tc)e 

- ln(l - Tc - p^{Tc)) + In {pi{Tc)e-P + P2{Tc)e-^^) 

- ln(l -Tc- P3{Tc) - P2{Tc)) + hi (pi(Tc)) - /? , 



if ;3 > -a/3 , 

if -a/2</3<-a/3 

if - a< 13 < -a/2 , 



(F6) 
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and 



(T) = , if /3 < -a 



(F7) 



i.e. the system is in non-invasive regime. 

In the limiting case of small synergy, |/3| ^ 1, the 
locus of critical points is given by a stright line, 



ln(l - T,) - {nb)l3 - - ln(l - 



n=l 

(F8) 

where (uh) is the mean number of bonds attached to 
the donor given at least one bond atttached. For square 
lattice with Tc = 1/2 and pi = P2 = 3/7, = 1/7, 
Eq. gives 



a ~ ln(2) - 12/3/7 



(F9) 



In the limiting case of strong interference, jS — 
— c», the phase boundary approaches the straight line 
asymptote, given by 

a = - ln(l - -p3(Te)) + In {pi{T,)e-^ + p2{T,)e"^P) 

P2{Tc) 



In- 



P3(re) 



--2/? = ln(6/5)-2/3 



(FIO) 



for square lattice. 

The locus of critical points for d-synergy given by 
Eqs (|F6| and (|F7| is shown in Fig. [16] (dot-dashed 
line). The intersections of the critical line with the 
straight lines, (3 = —a/n {n ~ 1,2,3), correspond 
to the changes in the regimes for the infection rates 
in Eqs. ()F6P and (|F7[) and thus lead to appearance 
of the kinks (discontinuities in the derivatives) on the 
critical line. Within our approximations, the values 
of pn for square lattice are such that there is only 
one kink on the phase boundary for simple analyti- 
cal model (corresponding to the crossing point of the 
phase boundary with (5 = —a/3) and in the asymptotic 
regime, /? — s- — cx), the phase boundary approaches lin- 
ear asymptote, /3 ~ /3q ~ a/2 (with positive constant 
P'o > 0). 

Kinks are also expected on the phase boundary for 
the exact model. In order to check this, we have 
analysed the behaviour of the exact phase boundary 
obtained numerically around the point of intersection 
with the line /3 = —a/2 (cf. Fig. [TBI) . We have tested 
the presence of a kink by fitting a linear dependence 
/3 = a + ba to the data above and below the intersec- 
tion. This procedure reveals a significant difference in 
the slope which takes values b = —0.97 and b = —1.26 
above and below the intersection, respectively. 

The locus of critical points in a — /3 parameter space 
derived within the simple analytical model (dot-dashed 
line in Fig. I16[) is similar in shape to that obtained nu- 
merically (dashed line marked by the circles in Fig. ll6p . 
However, due to simplifying approximations (i)-(iii), 
the analytical model does not capture the asymptotic 
behaviour, /3 ~ /3o — a (with positive constant f3o > 0), 
obtained numerically for /3 — >■ oo. This could be due 
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FIG. 16. Locus of critical points for d-synergy within simple 
analytical model (dot-dashed line) compared with the exact 
numerical data (dashed line marked by circles). The inter- 
sections of the critical lines with straight lines, f5 = —a/n, 
indicate the location of the kinks on the phase boundaries. 



to the fact that the actual probabilities p„ (possibly 
depending on a and /3 might be quite different from 
those used in the analytical model (see Eqs. (|F4|) ). The 
difference in gradients in the small-/3 limit between an- 
alytically and numerically found phase boundaries can 
be due to similar reasons. 



2. r-synergy 

For r-synergy, the transmission of infection from a 
donor to a recipient can occur in the presence of differ- 
ent number of infected neighbours (in addition to the 
donor) of the recipient (for concreteness, we consider 
a square lattice for which n = 1,2,3). Under assump- 
tions (i) and (ii), the mean transmissibility for the d-r 
pair is. 



(Fll) 



where g„ is the probability that the recipient has n in- 
fected neighbours (in addition to the donor) and syn- 
ergistic transmissibilities are given by Eq. (jF3p where 
n = 0,1,2,3. The probabilities g„ of different neigh- 
bourhoods of the recipient can be defined through the 
probability, t (parameter of the model), that a near- 
est neighbour of the recipient host (different from the 
donor) is in the infected state. 



C;li"(l 



(F12) 



where we used the assumption about the independence 
of infection events for different neighbours of recipient. 
Given assumption (iii), from Eq. (|Fip written in the 
form, 



Tl = 



(F13) 
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we can easily obtain the following equation for the locus 
of critical points: 



a 



ln(l - + In (qoit) + q,{t)e-P + q2{t)e-^^ + q^{t)e-^P) 
ln(l - - q^{t)) + In {q,,{t) + gi(<)e-^ + 92 ^e^'^) , 
ln(l - T, - 93 (t) - <Z2 W) + In (<Zo W + gi(i)e"'') , 

* ~ 9o(t)-Tc ' 



a (3> -a/3 

if -a/2<^<-a/3 

if -a<f3<~a/2 

if /? < -a . 



(F14) 



In the limit of large values of /? — ?■ 00, the value of a 
approaches the asymptotic value adoo), 

a ^ ae(^) = In . (F15) 

t -I c 

The value of adoo) can be both positive and negative. 
We know from exact numerical analysis that ac(oo) > 
{adoo) — 0.2) and thus we assume that 

go = (1 - <)' > 1 - , (F16) 

i.e i < 1 - (1 - Tc)i/3 (J- < ^j. square lattice). 
In the limiting case of small values of — > 0, 

3 

a ac(oo) = - ln(l -Tc) - Vn(7„/3 = 

1 l^^'J 

= - ln(l - T,) - , 

where (n) is the mean value of infected neighbours (ex- 
cluding the donor) of the recipient. 

In the limiting case, (3 —00, the behaviour of the 
analytical critical line depends on the value of t. In 
particular, if qo = {1 — t)^ > Tc then the critical line 
intersects the straight lines (3 — —a/n (resulting in 
kinks on the critical line) for all values of n = 1,2,3 
and becomes a vertical border at a = a* for (3 < —a*. 
Numerical data support such a scenario with the ver- 
tical border. The positions of both asymptotic value 
of adoo) and vertical border a, vary with the value of 
t. In Fig. [171 the analytical critical lines are shown for 
two values of the node occupation probability t. For 
t ~ 0.047, the behavior around small values of /? and 
position of the vertical border found numerically arc 
reproduced quite well by the analytical model. How- 
ever, the analytical curve strongly deviates from the 
numerical one for /3 3> 1 and fails to reproduce the 
value of adoo). If we try to mimic the value of adoo) 



by tuning t (see the curve for t = 0.152), then the gra- 
dient at small |/3| and position of the vertical border 
are significantly off the numerical values. Such devia- 
tions area consequence of approximations (i)-(iii) used 
in the analytical model. 

Overall, comparing the two types of synergy within 
the simple analytical model we can conclude that the 
main differences between them come from the presence 
of non-syncrgistic transmission events that are possible 
for r-synergy with any value of /3 when the recipient is 
challenged by a single infected neighbour. Such trans- 
mission events with transmissibility oc 1 — e^" are re- 
sponsible for the appearance of the asymptotic value 
adoo) and vertical border a^, for r-synergy. 
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FIG. 17. Locus of critical points for r-synergy within simple 
analytical model for two values of parameter t (as indicated 
in the legend) compared with the exact numerical data. 
The intersections of the critical lines with the straight lines 
P = —a/n indicate the locations of the kinks on the phase 
boundaries. 
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